Vortex states of rapidly rotating dilute Bose-Einstein condensates 
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We show that, in the Thomas-Fermi regime, the cores of vortices in rotating dilute Bose-Einstein 
condensates adjust in radius as the rotation velocity, fi, grows, thus precluding a phase transition 
associated with core overlap at high vortex density. In both a harmonic trap and a rotating hard- 
walled bucket, the core size approaches a limiting fraction of the intervortex spacing. At large 
rotation speeds, a system confined in a bucket develops, within Thomas-Fermi, a hole along the 
rotation axis, and eventually makes a transition to a giant vortex state with all the vorticity contained 
in the hole. 
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The energetically favored state of rotating superfluids 
such as^^He Uj and the alkali gas Bose-Einstein conden- 
is a triangular lattice of singly quantized 
Here we pose the question of the structure 
of the lattice and ultimate fate of the superfluid at ro- 
tation speeds so large that a sizable fraction of the fluid 
is filled by the vortex cores. A Type II superconduc- 
tor with a vortex array simply becomes normal at a high 
magnetic field, Hc2, where the vortex cores begin to over- 
lap 6]. While in a rotating Bose superfluid the cores of 
quantized vortex lines (if of fixed size) would completely 
fill the system at an upper critical rotation speed r2c2 
analogous to Hc2, a low temperature bulk bosonic sys- 
tem does not have a simple normal phase to which it can 
return, and thus the problem |7|. In He II the critical 
rotation rates at which the vortex cores approach each 
other are unobservably large, fJc2 ~ 10^^ rad/sec. In 
Bose condensed dilute atomic gases, by contrast, inertial 
forces can be comparable to interatomic forces, and the 
approach to tightly packed vortex lattices is within reach 
of experiment; indeed current experiments in harmoni- 
cally trapped gases achieve rotational velocities Vl that 
are a significant fraction of Oc2 |^ IE Hi IE 1^ - 

We consider a Bose gas at zero temperature in the 
Thomas-Fermi regime, with weak repulsive interactions 
described by an s-wave scattering length Os, and inte- 
grate out (in the sense of the renormalization group) the 
short range structure on scales of the vortex separation. 
The structure of the cores and the destiny of the vortex 
lattice at high rotational speeds depends fundamentally 
on the geometry of the container confining the fluid. We 
study here the modiflcation of the lattice induced by high 
rotational speeds in both harmonic traps and in cylin- 
drical hard- walled (square- well) buckets |3|- In a har- 
monic trap, where considerable work has been done on 
understanding the states of the fluid 0, 0| , the radial 
trapping potential, V±^ = ^muj\r'^ , where r is the per- 
pendicular distance from the axis of rotation, dominates 
the centrifugal potential. As the rotational frequency 
n approaches ajj_ , the system becomes quasi- two dimen- 
sional and should eventually enter a quantum Hall-like 



state 0, 0, Q, 0|; we do not treat this hmit here. 
We find that in a trap the cross-sectional area occupied 
by the vortex cores grows until they fill a limiting frac- 
tion, ~ 1/2, of the space of the system, and they never 
touch, even for rotational velocities arbitrarily close to 
This result goes beyond that of 0, in that we take 
into account (albeit approximately) the modification of 
the vortex core in the finite sized cell of the lattice, for 
given particle number in the cell. 

In a cylindrical bucket, on the other hand, where the 
centrifugal force tends to throw the fluid against the 
walls, we find that at a critical rotation speed, Sl/i = 
V2h/m(,oR, the fluid begins to develop a cylindrical hole 
in the center which grows with increasing Q; here R is 
the bucket radius, = h{2mgfi)~^^'^ is the coherence or 
healing length, g = Airfi^as/rn, and n is the uniform den- 
sity in a non-rotating cloud. In the bucket, the ratio of 
cross-sectional core area to the area of the unit cell of the 
(2D) vortex lattice is independent of for > and 
given by O^o/v^i?- A phase transition associated with 
core overlap at high vortex density is thus precluded for 
either trapping geometry. In addition, as we show, a sys- 
tem in a bucket eventually undergoes a phase transition 
in the limit of large 17 to a giant vortex phase with the 
vorticity concentrated along the cylinder axis. 

The scale of rotation velocities in a dilute gas in a 
container is set by fig = 2gn/h = h/m^Q ~ rtc2- For 
^0 ^ 0.2/Ltm nil, rio - 1-8 X 10^ rad/sec in ^^Rb and 
~ 6.9 X 10" rad/sec in ^SNa. 

We describe the ground state of the system with a 
rotating vortex lattice by an order parameter, ip{r), de- 
termined by minimizing the energy, E' = E — ^L^, in 
the rotating frame by means of a variational calculation; 
here Lz is the component of the angular momentum of 
the system along the rotation axis. With Ti = 1, 



E' = d 
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where V{r) is the trapping potential. For a single vor- 
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tex line of unit winding number on the cylindrical axis, 
i'{r) = f{r, z)e'^'^ in polar coordinates. To calculate E' 
we employ a Wigner-Seitz approximation replacing 
the triangular shaped cells of the (2D) lattice by cylin- 
drical cells of equal radius t, the assumption of equal 
cell area is consistent with observed lattices in traps 0. 
We assume straight vortex lines for simplicity, and that 
the lattice rotates at a uniform angular velocity re- 
lated to the (2D) density of vortices Uy — l/nP by 
0„ = Tmy/m = l/m£^; ^ly is in general smaller than 
the stirring velocity 57 by terms of order \/mE?. 

In general, ip in & cell is a unique function of the num- 
ber of particles therein, and can be constructed, e.g ., by 
the method of matched asymptotic expansion |l7| . It 
is sufficient here to use the simple approximation that 
within a cell, labelled by i. 



(a;/e)y^^(I), < a; < ^, 
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where x is the (polar) radial coordinate measured from 
the center r.; of the cell. To account for the density vari- 
ation in the system we allow the mean density n{ri,z) 
within a cell centered on r.i to depend on the radial po- 
sition of the cell and height. Ansatz ||2Jl is consistent 
with the exact boundary conditions, f{x — *■ 0,z) x, 
and df{x,z)/dx\x=i — 0. While the core radius in a cell 
should depend on the local density, for simplicity we take 
^ to depend only on the mean density in the system, an 
approximation adequate to bring out the essential physics 
of rotation at high speeds in the Thomas-Fermi regime. 
As one approaches the mean-field quantum Hall regime 
01 , it becomes necessary to include possible dependence 
of / on the trapping frequency as well. 
The mean density n{ri, z) in the cell is 



n(n 



dV|^P=n,(z)(l-C/2), (3) 



where C, = (C/-^)^ < 1 is the ratio of the core area to the 
transverse area of the unit cell of the lattice. 

Writing the local velocity asv = ^lyrd>+Sv to evaluate 
the kinetic energy in each cell, as in [l^], we find 



E = 



cPrn{r) [ ^-mnlr^ + {a + l)f7„ 



+V{r) + ^gn{r)b ) , (4) 

where we replace the sum over cells by a continuous in- 
tegration, with Ti ^ r. In Eq. the kinetic energy per 
cell is 7ra(^)n(r)/m and the interaction energy per cell is 
gn'^{r)b{C)n£'^/2, where 
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a(C) has a minimum at C — 0.56, while &(C) increases 
slowly with The first term in |0J is the energy of solid 
body rotation of the lattice, while the term ^ n{r)ily 
arises from (f) ■ 5v. The result Q is effectively 
the Thomas-Fermi approximation for the coarse grained 
structure, described by n{r), with the functions a and b 
accounting for the fine grained vortex structure. 

The total an gula r momentum in the Wigner-Seitz ap- 
proximation is 



d^rn{r) (l + Q.^mr'^) 



(6) 



the sum of the angular momentum per vortex cell and 
the solid body rotation contribution; is independent 
of a and b. In toto, 



E' 
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Minimization of Eq. lO with respect to fly implies that 
ily — — [a + l)N/I, where I = J d^rmr'^n{r) is the 
moment of inertia. For large systems, mR^Vl ^ 1, we 
neglect the small difference of fit, and fi; then 
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Minimization of Eq. (jS]) with respect to n{r) at fixed total 
particle number yields the Thomas- Fermi result. 



i{r) = [/i + - V{r)]/gb, 



(9) 



where is the chemical potential and fl = fi — afl. In a 
trap with transverse frequency lu±_ and longitudinal fre- 
quency LUz, the density has the usual form [l9j 



ntrap(r-) = (AM) {l~{r/Rtf 



with i?2 = {2il/m)/{Lj]_ - n'^) and Zf 
Then N = SwZtR'^fi/Wgb and 



i^^/Ztf) , (10) 
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By contrast, the density in a cylinder grows quadrati- 
cally with r, until at a critical velocity flh where /i = 0, 
the system develops a cylindrical hole in the center; in 
a trap the quadratic trapping potential dominates, pre- 
venting formation of a hole. Integrating Eq. Q over the 
volume of the cylinder, we find that the hole begins to 
develop at Q,, = V2bno^o/R = {2bD.o/mR^Y/'^ . Then 



ncyi(r) 



n + mri^ (r^ - i?2/2) /2gb, n < ft,,, 
mn^{r^ - Rl)/2gb, !7 > f)/,; 



(12) 
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where the radius of the hole is Rh = R\/l — The 
phase of the order parameter winds by 2TTmflRfi^ at the 
inner edge of the hole, corresponding to an array of phan- 
tom vortices of density n„ inside the hole. 

The energy per particle in the rotating frame is 
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The core size is determined by minimizing E' with re- 
spect to C at fixed N and O; differentiating Eq. © di- 
rectly we have 



da g{n)_db _ 
9C 2 dC~ 



(15) 



where (n) — Jn^ir)/ J n{r) is a mean of the density 
in the system. [Note that were we, more correctly, to 
solve for the vortex structure within each cell, then the 
integration in (n) in Eq. H15|) would be only over the given 
cell.] In a bucket the relevant regime is small where 
&(C) ~ l + C,/3, and da/dC, ~ — 1/2^, the minimum is thus 
at Ccyi = i^lg{n). For il < VLh, (n) = fi (l + \{n/VthY), 
and 



Ccyl - 

For Vl>Vlh, In) 



= ^nJl/il^, and 



fl < fl,, 
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(17) 



The total area occupied by the vortex cores thus scales 
as the ratio of zero temperature coherence length and 
system size and is independent of H, see Fig. QJi. 

In a trap, (n) = Afl/7gb; for small fl, <^ ^ 1 and 
da/dC, ~ — 1/2C, so that ^ grows with increasing il as 
{lo\ — Q?)^ , as in 0. However, as fl ^ uji_, ft and 
therefore (n) — > 0, and the solution of H15() requires 
da/dC, 0, which occurs at ^ ~ 1/2; this is the lim- 
iting value of Quite generally, the fact that (rt) 
as n — > uj± implies that the solution must minimize the 
kinetic energy in Eq. (Q. The cores can never overlap, 
and the system can, within mean-field theory, maintain 
a lattice of quantized vortices for all Jl < lu±. This be- 
havior is shown in Fig. ^p; in strongly elongated traps, 
as in a cylinder, the system reaches a self-similar regime 
in which the core size scales with the intervortex spacing. 

The rotation speeds needed to observe the increase of 
core size with n are in principle accessible within cur- 
rent harmonic trap experiments, prior to the onset of a 
quantum Hall state. In iV ~ 5 x 10'', lu±/uJz = 4.2, 



^/uj± < 0.7, and u}±/i}o = 0.008; an increase in aspect 
ratio by a factor 24, corresponding to the solid curve in 
Fig. should be sufficient to observe the strong in- 
crease of the core size for il < ujj_ . 
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FIG. 1: Vortex core area as a fraction of the Wigner-Seitz cell 
area: (a) ( vs. il/flo in a hard-walled bucket with R/^o = 500 
(solid) and R/^o ~ 1000 (dashed); (b) vs. fl/uj± in an 
anisotropic cylindrically symmetric harmonic trap with pa- 
rameters N = 5 X 10^, as/d± = 0.01, and lj±/uJz = 10^ (solid 
line), 10" (dashed), and 10* (dash-dotted). 



At sufhciently large il, a dilute gas in a cylindrical 
trap makes a transition to a giant vortex state, with 
vorticity concentrated in the center rather than spread 
throughout the system. Such behavior is seen in numer- 
ical simulations in traps that rise faster than harmonic 
[23,1^ (and also in the presence of vortex pinning (22l|). 
Dynamical formation of metastable giant vortices in har- 
monic traps with 7 to 60 units of vorticity is reported in 
Ref. [2^. At the level of Thomas- Fermi, which is ade- 
quate to describe the giant vortex regime in large systems 
{mR^n 1) a giant vortex, described by order param- 
eter, ip{r) = ■\/?iG(?')e""^; has an energy in the rotating 
frame, 



E' 



(fir nQ{r) 



2mr^ 



-nQ{r) — vil 



(18) 



Minimizing with respect to n{r) we find the analog of 
Eq. nQ{r) = {fl—v'^/2mr'^)/g, where /i = ^ + vVl. In 
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Thomas-Fermi, the giant vortex has a hole in the center 
of radius Rq = aR, where a = v / R^2mjl. In the Umit 
of large Q. ^h), where a ^ 1, we expand in powers 
of {VlQ/mR^n'^f/'^ and find 
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Comparing with the lattice energy, Eq. 1)14(1 . and using 
(|17|l . we find the transition to the giant vortex state at 



nr. = 



no 



no 



9(1 + 2a(9^/\/2i?)) 9(l + ln(i?/eo))' 



In this paper we have studied only the zero tempera- 
ture structure of rotating Bose gases in the long wave- 
length Thomas-Fermi regime. The full structure at fi- 
nite temperatures reflected in the phase diagram in the 
n-T plane (studied for a harmonically trapped gas in 
[2^ ') may be rich. For example, at sufficiently large rota- 
tion velocities and low temperatures such that thermally 
excited motion perpendicular to the boundary becomes 
frozen out, the fluid forced against the boundary of the 
bucket becomes effectively two-dimensional. Under these 
circumstances the fluid can possibly undergo a transition 
to a two-dimensional Kosterlitz-Thouless phase. 
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